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We show that hard-scattering factorization is violated in the production of high-pr hadrons in 
hadron-hadron colUsions, in the case that the hadrons are back-to-back, so that fcr factorization 
is to be used. The explicit counterexample that we construct is for the single-spin asymmetry 
with one beam transversely polarized. The Sivers function needed here has particular sensitivity 
to the Wilson lines in the parton densities. We use a greatly simplified model theory to make the 
breakdown of factorization easy to check explicitly. But the counterexample implies that standard 
arguments for factorization fail not just for the single-spin asymmetry but for the unpolarized cross 
section for back-to-back hadron production in QCD in hadron-hadron collisions. This is unlike 
corresponding cases in e'^e~ annihilation, Drell-Yan, and deeply inelastic scattering. Moreover, the 
result endangers factorization for more general hadroproduction processes. 
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I. INTRODUCTION 

The great importance of hard-scattering factorization 
in high-energy phenomenology hardly needs emphasis. 
Essential to its application and predictiveness is the uni- 
versality of parton densities (and fragmentation func- 
tions, etc) between different reactions. However, as can 
be seen from [l|, 0, [1, 01 , process-dependent Wilson lines 
appear to be needed in the inclusive production of two 
high-transverse-momentum particles in hadron-hadron 
collisions, i.e., in the process 



Hi + H2 — > i?.3 + -^4 + X. 



(1.1) 



In this paper we will show that this situation definitively 
leads to a breakdown of factorization. 

The standard expectation is that the cross section is 
a convolution of a hard scattering coefficient da, par- 
ton densities, fragmentation functions and a possible soft 
function: 



da 



d'^Pzd^p. 



^ j dcji+j^k+i+x ft/1 fj/2d3/kdi/i 
+ power-suppressed correction. 



(1.2) 

Here the sum and integral are over the flavors and mo- 
menta of the partons of the hard scattering, /j/// denotes 
a parton density, and du/i a fragmentation function. 

It is noteworthy that the classic published proofs for 
factorization in hadron-hadron scattering jH, Q only con- 
cerned the Drell-Yan process. There are a number of 
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difficult issues in the proof that are highly non-trivial 
to extend to other reactions in hadron-hadron collisions, 
even though Eq. (|1.2p is a standard expectation. 

We will examine the case that the produced hadrons 
are almost back-to-back. Then the appropriate factoriza- 
tion property is fcT-factorization, which entails 0] the use 
of transverse-momentum dependent (TMD) parton den- 
sities and fragmentation functions. However, the issues 
raised by our counterexample to factorization are suffi- 
ciently general that they create a need to examine very 
carefully the arguments for factorization in hadropro- 
duction of hadrons even in situations where ordinary 
coUinear factorization with integrated densities is appro- 
priate. In the case of fcT-factorization with TMD den- 
sities, the factorization formula needs the insertion of a 
soft factor S, not shown in Eq. (|1.2|) . 

The problems concern gluon exchanges between differ- 
ent kinds of coUinear line, as in Fig. [7] below. To obtain 
factorization, the gluon attachments must be converted 
to Wilson lines in gauge-invariant definitions of the par- 
ton densities and fragmentation functions. This relies Q 
on the use of Ward identities applied to approximations 
to the amplitudes. But the approximations are only valid 
after certain contour deformations on the loop momenta. 

Bacchetta, Bomhof, Mulders and Pijlman [l|, [3, H, 0| 
argued that because of the complicated combination of 
initial- and final-state interactions, the Wilson lines must 
be modified. What is not so clear is the interpretation of 
their result. So in the present paper we present an argu- 
ment to make fully explicit the failure of factorization. 

Since the issue is one of factorization in general, and 
not just specifically in QCD, we clarify the issue by ex- 
amining a particular process in a model field theory. The 
process is a transverse single-spin asymmetry of the kind 
controlled by a Sivers function. This is a case where prob- 
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lems in the contour deformation directly affect the value 
of the cross section at the lowest possible order of pertur- 
bation theory. Our model field theory is simple enough 
that the calculations and their interpretation as imply- 
ing factorization violation arc unambiguous. But, as we 
will explain in the final section, we expect the failure of 
factorization to be more general: our particular process 
and model simply make it very easy to see the failure. 

II. CONSTRUCTION OF MODEL 

Since proofs of factorization apply to quantum field 
theories in general (if they are renormalizable) , the con- 
struction of a counterexample, to demonstrate and to 
understand a failure of the normal methods of proof, is 
conveniently done in a simple model theory. 

Our model resembles the one used by Brodsky, Hwang, 
and Schmidt [8] in their discussion of single spin asym- 
metries. It is defined as follows: 

• The gauge group is abelian. This simplifies the 
graphs, and allows the next feature. 

• The gluon is massive. This avoids the discussion 
being confused by actual infra-red divergences in 
the S-matrix. 

• The initial-state particles correspond to Dirac fields 
that are neutral under the gauge group. We will 
call them hadrons. The fields need to be Dirac 
fields in order to have the single transverse spin 
asymmetries that we will examine. We will use two 
types of hadron. 

• Each "hadron field" Hi will have a coupling to a 
Dirac field ipi and a scalar field (f>i, which in 
would be called a diquark field: 

Xi(^H,^P,(bl+^,H,(b?j . (2.1) 

The quark field ^pi has a coupling gi to the gauge 
field, and the scalar field (pi has the opposite cou- 
pling. 

• All the masses in the theory are comparable, to 
avoid confusing the calculation with logarithmic de- 
pendence on large ratios of masses. 

In our analog of hadroproduction. Fig. [71 the two initial 
state hadrons are of the two different types, and we use 
this to simplify our argument for non-factorization. The 
lines in the the lower part of graphs are chosen to be 
those of gauge coupling gi, while the lines in the upper 
part of the graph are those with coupling 52- But if 
the attachments of the gluon to the upper lines were in 
some way to correspond to a Wilson line in the parton 
density in the hadron in the lower part of the graph, the 
charge would have to be gi . Since gi and g2 are arbitrary, 
there is no way to make a correspondence between the 



graph and the Wilson line formalism, provided that the 
contribution is non-zero, as we will demonstrate. This 
will also insulate us against sign errors and the like. 

We will also choose the detected outgoing particles H3 
and H4 to correspond to the scalar fields. The sole pur- 
pose here is to simplify the Dirac algebra slightly, thereby 
making the calculations more transparent and elemen- 
tary. 

One feature of our counterexample appears to be very 
special to an abelian gauge theory. This is that the two 
couplings gi and 32 need have no relation to each other: 
there is a continuous infinity of representations of the 
gauge group. In contrast, there is a single value of the 
coupling g for all the fields in a non-abelian theory. The 
role of the ratio 52 / gi is now taken over by the represen- 
tation matrices for the different fields in QCD (triplet, 
antitriplet, octet), with the different couplings related 
by factors of rational numbers. So in any particular ex- 
ample there is a potential for a numerical coincidence 
between the sizes of the numerical values of the graphs, 
which could then appear to give consistency with factor- 
ization. A counterexample to factorization would then 
be more complicated, with a comparison of cases with 
different kinds of partons (quarks, antiquarks or gluons), 
cf. 

III. REVIEW OF SIDIS AND DY 

We now review how 0] a transverse single-spin asym- 
metry (SSA) arises in semi-inclusive deep-inelastic scat- 
tering (SIDIS), at the level of onc-gluon-exchange, and 
how it determines the Wilson line that defines parton 
densities. Then we review the differences that give factor- 
ization with an exact sign reversal in the Sivers function 
for the Drell-Yan process 0, [l^ . This will give us meth- 
ods of calculation that will give us a very elementary way 
to obtain the SSA for the process (jl.ip . 

A. SIDIS 

With the electromagnetic part of the scattering fac- 
tored out, SIDIS is the process 

riq)+Hip)^H'{r)+X. (3.1) 

We use light-front coordinates in which the incoming mo- 
menta are 

(3,2) 

The detected outgoing particle is defined by a longitu- 
dinal momentum fraction z and a transverse momentum 
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We will assume that Q is large and that the detected The lowest-order graph Fig. [T] gives the following con- 
transverse momentum rj- is of order a hadronic mass tribution to the differential structure tensor: 
scale m. 



dW^'•' 
dz d^rT 



A2 
47r 



dk+dk- 
~J2^ 



k-+q-\ (2fc^ + g'')(2fc'' + ?'•') 



{2nf5{{q + kf ~ ml) 5{{p - kf - ) x 



1. 



X -Tr[(|* + mn) (1 + 75^) {f - ^ + m^)] 



(3.4) 



The internal partons are all coUinear to the target, i.e., k^ ^ p'^ , k^ ^ rn^/p^, kr ~ to, and to leading power in 
Q, parton model kinematics apply, so that k^ ~ xp'^ and z ~ 1. We will assume throughout that the spin vector 
s corresponds to a transverse spin (in the {QtP) frame), and that it is normalized so that its extremal value obeys 

= —1. Since there is only one initial-state hadron, we do not bother with labels to indicate the kind of hadron 
(e.g., Ai or A2). 

The above formula is simply related to a parton density: 



(p^ - q^'p ■ q/q^){p'' - q^p ■ q/q^) 



dz d^TT p ■ q 

where the parton density at lowest order is 



Pq{x, kj 



5{z — 1)P^/h{x, tt) + power-suppressed correction, (3.5) 
iTr[(|* -I- niH) (1 + 75;^) {f - ^ + m^)] 



167r3 



(3.6) 



FIG. 1: Lowest-order graph for SIDIS in our model. The ini- 
tial state particle is a color-singlet Dirac particle. The spec- 
tator line is for a Dirac "quark" field, and the active parton is 
for a scalar "diquark" field, as can be enforced by a condition 
on the detected outgoing particle of momentum q + k. The 
arrows on the internal lines indicate the flow of color charge, 
and the vertical line cutting the graph denotes the final state. 



as follows from the conventional operator definition. 

There is in fact no polarization dependence in this or- 
der, i.e., the Sivers function vanishes. At a simple calcu- 
lational level, this occurs for two reasons. One is that 75 
gives a non-zero contribution only when multiplied by at 
least four regular Dirac matrices. 



[fcy + "1^(1 — x) + m'^x — ■m'jjx{l — x)] ' 
I 



reason for the vanishing is that the trace p.7p is imag- 
inary, while the rest of the graph is real, so that the 
contribution to a cross section must be zero. 

The lowest-order graphs for a non-zero SSA are the 
one-gluon-exchange graphs in Fig. ^ which get an imagi- 
nary part from an intermediate state that can go on-shell; 
this state is made by the lines with momenta q + k — I 
and p — k + I. Standard power-counting shows that the 
exchanged gluon can only be coUinear to the target or 
soft. The minus momentum of the gluon is trapped in 
the region l~ ~ ni? /p'^ by the other target-coUinear lines. 
The on-shell intermediate state corresponds to small an- 
gle elastic scattering, and so to very small Z+, of order 
p^w? IQ"^. There the only significant dependence on 
is in the upper parton propagator. Multiplied by the 
neighboring gluon vertex this gives 



-g{2q^' + 2fc'^ - /^) 



tr75 



= A^e,^^,,a^h^c^'d'', 



(3.7) 



while in (|3.6p there are at most three. This reason for 
a vanishing SSA would no longer apply in a more com- 
plicated model or with higher order graphs. The second 



{q + k — l)^ — -I- ie —1+ + other terms -I- le' 

(3.8) 

where the "other terms" are small or independent of Z+, 
and we have taken a leading power approximation for the 
momenta in the numerator. The contour of integration 
of Z+ can therefore be deformed into the lower half plane 
until I is target-collinear. In that case the "other terms" 
in (|3.8[) are negligible, and the denominator can be re- 
placed by its eikonal approximation l/(— + ie). This, 
together with a leading-power approximation in the nu- 
merator, shows 9] that the gluon exchange correction is 
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FIG. 2: Virtual one-gluon-exchange corrections to Fig. [T]that give a SSA. 
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FIG. 3: Virtual one-gluon-exchange correction to parton den- 
sity. The upper double line is the Wilson line, and the graph 
shown, together with its hermitian conjugate gives the first 
contribution to the Sivers function. 
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FIG. 4: Lowest-order graph for the Drell-Yan process. 



equivalent to a contribution to the parton density with a 
suitable Wilson line, Fig. [31 

Let us perform the and k~ integrals by the on-shell 
conditions on the final state, and let us perform the l~ 
integration by contour integration. Then the necessary 
imaginary part comes simply from the imaginary part of 
p.Sp and thus from the replacement 

This gives rise to an SSA with the aid of the trace 

iTr[(|i + 771^)75^ {f-$ + t + m^)l^ {f - $ + m^)] 
~ 2iejksn''p+[mH{l -x)+ m^,], (3.10) 

where the approximation, good to leading power, arises 
from the neglect of the small components of I with respect 
to the transverse components. The two-dimensional e 
tensor obeys ei2 = 1. Since the denominator in the in- 
tegrand is not azimuthally symmetric in It, the integral 
over / gives a non-zero result for the SSA from the whole 
graph. 

The two graphs in Fig. [5] are related by hermitian con- 
jugation and so they give equal contributions to the SSA. 

B. Drell-Yan 

The Drell-Yan (DY) process, 

HA{s)+HB^l*{q)+X, (3.11) 



is treated quite similarly. We examine the cross section 
differential in qT, and investigate a possible SSA, with 
Ha having a transverse spin vector s. 

In our model the lowest-order graph is Fig. 3) It is 
readily shown to be the (convolution) product of two 
transverse-momentum-dependent parton densities, and, 
just like the SIDIS process, it has no SSA at this order. 
Note that to have the process occur at the order shown 
within our model, the initial-state hadrons Ha and Hb 
must be antiparticles of each other. 

The imaginary part in the amplitude needed to get an 
SSA with respect to the transverse spin s of the lower 
hadron Ha arises from graphs such as those in Fig. [S] 
Graphs (a) and (b) work just like those in Fig. [2] for 
SIDIS, except that the gluon couples to the incoming 
scalar antiparton instead an outgoing scalar parton, with 
the necessary reversal in sign of the coupling. Thus Eq. 
(|3.8[) is replaced by 

(fcs + l)-^ — TO^ -\- ie ^+ — other terms -I- it 

which gives an imaginary part exactly opposite to that for 
SIDIS. The relative sign of the Z+ term and the ie is now 
that for an initial-state interaction, so that the Wilson 
line in the operator definition of the parton density is now 
past-pointing instead of future-pointing [9;]. As shown 
in |9|, an exact reversal of sign of the Sivers function 
between SIDIS and DY follows from the time-reversal 
symmetry of QCD. 

No contribution to the SSA is given by graphs, like Fig. 
[SJc), in which the both ends gluons couple the the active 
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FIG. 5: Virtual one-gluon-exchange corrections to Fig. |4]relevant for a SSA when the lower hadron has transverse spin s. Graph 
(a) and its hermitian conjugate (b) have imaginary parts (at the amplitude level) that give a non-zero SSA. Gluon exchange 
between the active partons, graph (c) and its not-shown conjugate, gives an imaginary part in the vertex correction that does 
not give a SSA. Spectator-spectator gluon exchange graphs, (d) and (e), do contribute individually to the SSA, but the two 
contributions cancel (at leading power). 



partons, or where the gluon couples the upper spectator 
quark to the lower active parton. In our model this is 
trivial: there are too few Dirac matrices on the lower 
line to give spin-dependence. In a more general case, the 
annihilating partons could be Dirac fields, and then the 
lack of spin-dependence arises because of the eikonaliza- 
tion of those parton lines on the lower side of the graph 
that would contribute to the imaginary part. 

But spectator interaction graphs, (d) and (e), do indi- 
vidually contribute to the SSA. The imaginary part arises 
from putting all four spectator lines on-shell. Provided 
the gluon momentum is routed the same way in both 
graphs before integration, e.g., to the left as shown, the 
two contributions are equal, but exactly opposite in sign. 
The cancellation is exactly one of those needed to prove 
factorization Q , and involves a sum over all allowed cuts 
of a particular graph. 

However, the two graphs have final states with parti- 
cles of different momenta. This accounts for the differ- 
ence compared with the non-cancellation of the graphs 
in Fig. [2] for SIDIS, where there is also a sum over cuts. 
For the SIDIS graphs there is a requirement on the trans- 
verse momentum of the struck parton. In contrast, the 
cancellation between the two spectator-spectator graphs 
in the DY process occurs because our cross section is 
the fully inclusive DY cross section; no requirement was 
placed on the final-state in the target fragmentation re- 
gions. (As is well-known from the case of diffractive hard 



scattering, both theoretically [Tll,[T2| and experimentally 
[13I . [Til , factorization fails when target-relative restric- 
tions are imposed.) 



IV. HADROPRODUCTION OF HADRONS 

We now have the tools to make an extremely stream- 
lined construction of a counterexample to factorization 
for the process of hadro-production of high transverse 
momentum hadrons. Hi + H2 — > H3 + H4 + X. We 
again use an SSA because non-factorization occurs with 
one gluon exchanged beyond the lowest order in which 
the reaction occurs at all. We choose Hi to be the po- 
larized hadron, and we choose the hadrons Hi and H2 
to be of the two different flavors in our model. We also 
choose the detected final-state particles H3 and H4 to 
correspond to the two flavors of scalar parton. The high- 
transverse-momentum particles H3 and H4 are chosen 
to be almost back-to-back azimuthally (relative to the 
collision axis), so that transverse- momentum-dependent 
parton densities and fragmentation functions are needed 
for describing a factorized cross section. 

The single lowest-order graph for the process is shown 
in Fig. [6] Its hard scattering is just the gluon-exchange 
subgraph. The cross section is the convolution of the 
hard scattering with a transverse-momentum-dependent 
parton density in each hadron. The fragmentation func- 
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FIG. 6: Lowest order graph in model for hadroproduction 
of hadrons of high transverse momentum. The initial state 
particles are color-singlet Dirac particles. The spectator lines 
are for Dirac "quark" fields of charges g\ and g2, and the 
active partons are for scalar "diquark" fields. The exchanged 
gluon line is thickened to denote the hard scattering. 



tions in this order are trivial delta functions. Although 
the longitudinal momenta of the incoming partons for 
the hard scattering are determined from the kinematics 
of i?3 and i?4, only a sum of their transverse momenta 
is determined. Hence a convolution over the transverse- 
momentum densities is needed. As before, there is no 
SSA at this order. 

The graphs giving the lowest-order SSA are shown in 
Fig. [T] They have an extra gluon exchanged between 
the spectator line in the polarized hadron and one of the 
active partons. As with the DY process, a sum over cuts 
of graphs with a spectator-spectator interaction cancels, 
while exchanges purely between active partons give no 
SSA. An exchange purely between the spectator {pi —ki) 
and the active parton line (/ci) in the polarized hadron 
has no relevant intermediate on-shell state and therefore 
does not contribute to the SSA. 

So only the three graphs shown in Fig. [7| contribute to 
the SSA, together with their conjugate graphs. Exactly 
as in our discussion of SIDIS and DY, the graphs are the 
same to the leading power except for eikonal factors from 
the parton lines connecting the upper end of the gluon 
to the hard scattering. The appropriate replacements for 
these lines are 



-.g2(2fc^ - 
(fc4 — ly — + 

(fc2 + 0^ - "^02 + 
-ffl(2fc3^-?^) 



(4.1) 
(4.2) 
(4.3) 



for a total of iTT{gi + 2g2)S{l'^). 

The (?i term corresponds to a gluon coupling to a 
future-pointing Wilson line in the operator definition of 
the parton density, the same as for SIDIS. However, it 
is impossible for the contribution proportional to g2 to 



be represented in terms of a Wilson line connecting the 
two parton fields for the distribution of partons of type 
(pi in the hadron Hi. This is simply because the cou- 
pling for any such Wilson line has to correspond to the 
color charge of the parton, i.e., gi, and not g2. The full 
Wilson line, or some generalization, is needed because ex- 
changes of multiple gluons also contribute. The quantity 
we are looking at is definitely associated with the hadron 
Hi rather than being in some kind of exotic soft factor, 
since the attachment of lower end of the gluon line to the 
spectator parton is necessary for the non-zero SSA; the 
triviality of this observation is a special feature of our 
particular model. 

The contribution to the SSA is therefore non-universal 
and does not correspond to a parton density. That is, 
factorization is broken. 

Of course, the fact that the contribution is obtained 
from an eikonalized line indicates that it can be obtained 
from some kind of representation in terms of Wilson line 
operators. But the matrix element is for some more com- 
plicated and non-universal object P, 0, H, 0] that cannot 
be treated as a par_ton density. It is allied to the objects 
used by Balitsky [l5| to discuss scattering at high-energy 
and small angles. The eikonalization indicates that sub- 
stantial simplifications are possible. But that situation 
would go well beyond normal factorization. 



V. DISCUSSION 

We should first emphasize that there is a large over- 
lap between the present paper and the work in Refs. 
[Il 0, S 01 • What is not so clear from the earlier work 
is whether factorization in any standard sense continues 
to hold for in the process (|l.ip . For example, in [l|, 
we read "We have assumed factorization to hold in this 
treatment of TMD effects although it is, at present, cer- 
tainly not clear whether such a factorization holds for 
hadron-hadron scattering processes with explicitly TMD 
correlators." 

Our primary result is to show by a counterexample 
that hard scattering fcr-factorization with universal par- 
ton densities fails for the production of high px hadrons 
in hadron-hadron collisions, when a pair of measured 
hadrons is close to back-to-back azimuthally. The over- 
all issue is that in a gauge theory arbitrary exchanges 
of gauge fields between different coUinear groups ("jets") 
can occur without any power suppression. To obtain fac- 
torization it is necessary to show that the sum over these 
exchanges can be absorbed into the definitions of the par- 
ton densities and fragmentation functions, assisted by 
certain cancellations. A full proof will be quite general, 
applying to a general gauge theory and to many reac- 
tions. So one particular counterexample is sufficient to 
show that such a proof does not exist; we can then choose 
the counterexample for maximum clarity and simplicity. 

Even for those cases where factorization does hold, the 
need to make suitable definitions of the parton densities. 
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(a) (b) (c) 

FIG. 7: One-gluon exchange in model for hadroproduction of hadrons of high transverse momentum. The specification of the 
process is as in Fig. |6] Only graphs that contribute to the SSA are shown. Hermitian conjugates of these graphs also contribute, 
with an equal value. 



etc, SO as to absorb the effects of the gluon exchanges 
indicates that the parton densities, etc, can always be 
regarded as effective densities [l6| . The primary prac- 
tical issue is whether they are universal, i.e., the same 
for all reactions. In a certain sense, the well-known scale 
dependence of the densities is a kind of non-universality: 
different parton densities are needed when the scale of 
the hard scattering is given a large increment. But there 
is an evolution equation for the scale-dependence, and 
this applies to an individual parton density. No details 
or specification of the hard scattering is needed to treat 
the evolution equation, either to derive it or to apply it. 
We should therefore refer in this case to "modified uni- 
versality" , not to non-universality. Similarly the reversal 
of the sign of the Sivers function between SIDIS and DY 
processes is a case of modified universality. 

At the upper end of the exchanged gluon in our coun- 
terexample, the interactions can be treated in the eikonal 
approximation. This is very similar to other discussions 
of partons passing through the gluon field of another 
hadron. A selection of relevant papers is 
Much of that work concerns the small x region, diffractive 
scattering, etc, whereas our counterexample applies in 
the fully conventional region where normal parton-model 
concepts are generally considered as fully applicable, i.e., 
parton fractional momenta are moderate and the scale of 
the hard scattering is comparable to the center-of-mass 
energy rather than being much less. 

Of course, interesting simplifications do occur, so that 
useful quantitative estimates can surely be obtained for 
the non-factorizing effects. However the methods are 
rather different than those for conventional factorization. 
Refs. [H, [13, H [ill indicate that the effects of the eikon- 
alized interactions are substantial, so that the numerical 
effects of non- factorization should be significant; in the 
present paper we did not estimate the numerical size of 
the non-factorization. 

The gluon exchanges in our counterexample are clearly 
tied to the target hadron at their lower end. But the 
coupling at the upper end concerns some parton other 



than the one coming out of the lower hadron. The non- 
canceling terms are sufficiently tied to the color flow at 
the hard interaction that they are not universal in any 
normal sense. This is the clearest indication of non- 
universality. 

The reader should not be misled by specific features of 
our counterexample into supposing that the failure of fac- 
torization is correspondingly restricted. These features 
include: the use of an SSA, the particular features of the 
model, and the particular order of perturbation theory. 
The use of the SSA is simply a way of getting the maxi- 
mal conceptual sensitivity to problems in constructing a 
proof of factorization. For an unpolarized cross section, 
we would need an extra gluon to be exchanged in order for 
the nonfactorization issues to arise, from graphs such as 
those in Fig. [S] Evidently, to demonstrate nonfactoriza- 
tion explicitly in this case, the number of graphs would be 
larger than in our example, and the explicit calculations 
would be much more lengthy. Standard power-counting 
arguments show that the contribution of this and related 
graphs is of leading power. It is very important to de- 
termine whether or not the sum of the potentially non- 
factorizing contributions actually does or does not cancel 
in the unpolarized cross section. 

Similarly, the choice of quantum numbers of the parton 
fields, of the abelian gauge group, and of the quantum 
numbers for the detected particles was simply to provide 
maximum transparency and simplicity to the counterex- 
ample. 

The fact that non-factorization can only occur with 
at least two extra ghions in a unpolarized cross section 
might suggest that the non-factorization is at high order 
in the strong coupling and therefore substantially sup- 
pressed. However the region of interest is at low virtual- 
ity for the gluons, so that the appropriate coupling is for 
a low momentum scale, where QCD is a strongly coupled 
theory. 

Even so, the number of extra gluons needed implies 
that the effects of non-factorization will only appear in 
quite high order in conventional perturbative QCD cal- 
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FIG. 8: The exchange of two extra gluons, as in this graph, 
will tend to give non-factorization in unpolarized cross sec- 
tions. 




FIG. 9: In a conventional perturbative QCD calculation for 
an unpolarized partonic cross section, non-factorization by 
the mechanisms discussed in this paper would first appear in 
graphs of this order. 

culations. Normally one performs calculations with on- 
shell massless quarks and gluons, and extracts coUinear 
divergences that are grouped with parton densities and 
fragmentation functions; any remaining divergences can- 
cel between graphs. Non-factorization in the hadronic 
cross section corresponds to uncanceled divergences in 
quark-gluon calculations. The lowest order in which the 
mechanisms we have discussed could possible give an 
uncanceled divergence in unpolarized partonic cross sec- 
tions is next-to-next-to-next-to- leading-order (NNNLO), 
as in Fig. [HI The region for the uncanceled divergence is 
where the lower gluon is coUinear to the lower incoming 
quark, and two of the exchanged gluons are soft. This 
graph is at least one order beyond all standard pertur- 
bative QCD calculations. 

Because our calculations directly concern cross sec- 
tions that use transverse-momentum-dependent parton 
densities, a certain amount of care is needed in inter- 
preting the results. The natural direction for the Wilson 
lines is light-like, as from Eq. (|3.8p . However light-like 
Wilson lines give divergences in transverse-momentum- 
dependent densities These are due to rapidity di- 
vergences po| in integrals over gluon momentum; they 
cancel Q in conventional parton densities only because 



of an integral over all transverse momentum in integrated 
parton densities. The solution adopted by Collins, Soper 
and Sterman 0] (CSS) was to define parton densities 
without Wilson lines but in a non-light-like axial gauge. 
The gauge-fixing vector introduces a cut-off on gluon ra- 
pidity, and then an evolution equation with respect to 
the cut-off was derived. The non-perturbative functions 
involved in this CSS evolution equation have been mea- 
sured (e.g., [l^l) in fits to DY cross sections, and would 
be an essential ingredient in testing non-factorization. 

However, there are some unsatisfactory features of the 
use of axial gauges, which are made particularly evident 
in polarized cross sections. This includes complications 
concerning gauge links at infinity [2^ , when a Wilson line 
formalism is used. A much better definition is to use a 
non-light-like Wilson line. This again obeys an equation 
of the CSS form. It is also possible to use a subtractive 
formalism [2^ [2^ with light-like Wilson lines but with 
generalized renormalization factors involving vacuum ex- 
pectation values of Wilson lines, which also implement a 
rapidity cutoff, and lead to a CSS equation. 

To test the predicted non-factorization, we simply need 
predictions of high-pr hadrons in hadron-hadron colli- 
sions, made on the basis of fits to parton densities in 
DIS and DY and to fragmentation functions in e"'"e~ and 
SIDIS 24]. Probing the SSA would be particularly in- 
teresting, and such measurements are underway at Rel- 
ativistic Heavy Ion Collider (RHIC) [H, [H. The same 
physics is probed in the transverse shape of jets, and 
would be worth investigating. 

Our counterexample applies in a kinematic region 
where the normal intuitive ideas of the parton model 
appear quite appropriate, even with a generalization to 
/cT-factorization. Therefore it forces us to question un- 
der what conditions factorization is actually valid and to 
what extent it has actually been demonstrated. It cannot 
be assumed that naive extensions of apparently estab- 
lished results are applicable beyond the cases to which 
the actual proofs explicitly apply. 

For hadron-hadron collisions, factorization has been 
proved [sl, ^ for the Drell-Yan process integrated over 
transverse momentum or at large transverse momentum 
(of order Q). These proofs apply in the presence of gluon 
exchanges of the kind that we discuss in the present pa- 
per. But these papers do not go beyond this, to the pro- 
duction of hadrons. Because factorization is important to 
all aspects of hadron-coUider phenomenology, it is critical 
to solve this problem for the hadroproduction of high-p^ 
hadrons. Given our counterexample to fcT-factorization, 
a proof of factorization can only succeed in a situation 
where conventional coUinear factorization is appropriate. 
For dihadron production this is when the hadron-pair has 
itself large transverse momentum or when the pair's out- 
of-plane transverse momentum is integrated over a wide 
range. 

In fact, Nayak, Qiu and Sterman 27] have recently 
given strong arguments that coUinear factorization does 
indeed hold in such a situations. The graphs examined 
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are similar to ours. They apply Ward identities to prove 
an eikonalization generalizing our specific calculations. 
Then they observe that a unitarity cancellation occurs 
of a kind endemic in factorization proofs [1, [2^. This 
concerns graphs that are related by different placements 
of the final-state cut. In our model, one example is given 
by Fig-IHa) and (b), and another is Fig. [71^ a) and its con- 
jugate. Such cancellations fail in our examples, because 
the final-states of the related graphs have different trans- 
verse momentum, and the cross section is not sufficiently 
inclusive in transverse momentum. 

Mechanisms that cause fc-r-factorization to fail in back- 
to-back hadron production also tend to cause resumma- 
tion methods to fail. They will also tend to break factor- 
ization or cause large perturbative corrections when de- 
tailed distributions of final-state hadrons are examined. 
Since many such cases are implicit in the analysis of com- 
plicated multi-jet cross sections, and of jet shapes and 
the like, encountered in searches for new physics, further 
understanding is essential as are quantitative estimates 
of the effects. They can have a particularly important 
effect in the extrapolation to the LHC of quantitatively 
measured distributions at the Tevatron, for example as 
embodied in Monte-Carlo event generators. The methods 



of [IB, ^^1^ be important. Probably some of these 

effects have already been modeled in some approximation 
and in at least some Monte-Carlo event generators, for 
example by the soft color interaction model [29j . 

Troublesome though it may be for phenomenology, 
breaking of factorization should be viewed not as some 
kind of failure, but as an opportunity. Examination of 
the distribution of high-transverse-momentum hadrons 
in hadron-hadron collisions will lead to interesting non- 
trivial phenomena. 
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